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ABSTRACT 


This  report  documents  the  mathematical  model  of  a  combined 
maneuvering,  stationkeeping,  and  seakeeping  simulation  computer  program  for 
two  different  classes  of  mine  hunting  vessels.  The  model  features  a  mix  of 
physical  based  models  for  the  propulsion  systems  and  control  surfaces  and 
hydrodynamic  coefficient  based  hull  model.  The  approach  follows  the  concept  of 
the  modular  maneuvering  model.  A  quasUsteady  assumption  is  utilized  for  the 
seakeeping  motion  effects.  This  allows  the  calm  water  maneuvering  motions  to 
be  calculated  separately,  and  the  six  degrees  of  freedom  linear  responses  in 
waves  to  be  superimposed.  The  effects  of  wind,  current,  and  second  order  mean 
drift  forces  are  included. 

ADMINISTRATIVE  INFORMATION 

This  work  was  funded  by  the  Coastal  Systems  Station,  Panama  City,  Florida,  under  Work  Request 
WR51 1991  and  is  identified  at  CDNSWC  by  Job  Order  Number  1-1561-432. 

INTRODUCTION 

The  Carderock  Division  of  the  Naval  Surface  Warfare  Center  (CDNSWC)  was  originally  tasked  in  1992  by 
the  Coastal  Systems  Station  (CSS),  Dahlgren  Division,  Naval  Surface  Warfare  Center,  in  Panama  City,  Florida  to 
develop  a  simulator  model  for  ship  maneuvering  and  stationkeeping,  including  seakeeping  effects.  That  simulator 
was  to  be  incorporated  into  a  team  training  simulator  for  the  MCM-1  (Avenger)  class  of  mine  hunter  vessels. 
Training  scenarios  would  cover  all  aspects  of  the  MCM  mission;  mine  sweeping,  neutralization,  and  transiting, 
possibly  in  the  presence  of  wind,  waves,  and/or  a  current.  The  result  of  this  effort  was  called  Version  1  of  the 
MCMSIM  computer  program,  as  documented  in  [1], 

In  1995,  CSS  tasked  CDNSWC  to  expand  the  original  simulator  model  to  allow  simulations  of  the  MHC-51 
(Osprey)  class  of  mine  hunters.  The  primary  modification  was  to  be  the  change  of  the  propulsion  system  from 
controllable  reversible  pitch  propellers  to  twin  vertical  axis  propellers.  Other  changes  and  modifications  include  the 
addition  of  roll  to  the  maneuvering  equations  of  motion,  a  physically  based  engine  model,  real  time  graphical  display 
of  simulation  results,  and  the  option  of  automatic  execution  of  standard  maneuvers  (turning  circle,  spiral,  and  zig¬ 
zag).  The  result  of  the  modifications  to  the  computer  program  (and  underlying  model)  is  now  called  MANSIM 
Version  2.  This  document  is  meant  to  replace  reference  [  1]  in  its  entirety. 

Originally,  the  basic  physical  requirements  for  the  combined  maneuvering  simulators  were  to  operate  in 
real  time  on  an  Intel  80386/80387  personal  computer  operating  at  25  MHz,  and  be  written  in  the  FORTRAN 
programming  language  using  no  machine  specific  modules.  For  the  development  of  MANSIM  Version  2,  an  Intel 
80486  DX  50  based  computer  was  used  with  a  32bit  FORTRAN  compiler  running  on  the  Windows  NT  operating 
system.  On  this  system,  a  minimum  of  five  to  ten  times  real  time  performance  has  been  obtained  (depending  on 


program  execution  options).  The  eventual  platform  for  the  simulator  program  is  a  Silicon  Graphics  model  4440 
workstation  operating  at  40MHz.  The  source  code  will  be  translated  by  CSS  into  the  ADA  programming  language 
for  final  inclusion  into  the  team  trainer  system.  In  the  final  form,  the  simulation  will  be  running  at  a  16Hz  sample  rate 
with  approximately  20  to  30  milliseconds  allowed  for  each  time  step. 

For  the  original  work,  two  different  starting  points  were  evaluated.  First,  one  of  the  existing  maneuvering 
programs  could  have  been  used  and  extensively  modified  to  support  the  zero  speed  special  case  required  for 
stationkeeping  [2,  3,  4].  Second,  a  new  program  could  have  been  written  from  scratch.  Either  option  needed  to  be 
coupled  with  a  wave  induced  ship  motion  calculation  methodology.  It  was  determined  that  because  of  the 
requirement  of  intrinsic  support  for  zero  speed  as  well  as  backing  and  sidewise  motions,  a  new  program  should  be 
written.  This  also  allowed  the  new  program  to  be  written  in  a  streamlined,  minimal  form  for  maximum  computational 
speed. 

The  development  of  the  original  simulator  model  represented  the  first  opportunity  at  CDNSWC  to  develop 
a  maneuvering  model  based  on  the  modular  maneuvering  model  concept.  The  modular  model  concept  differs  from 
the  traditional  approach  in  several  ways.  The  traditional  approach  is  to  use  hydrodynamic  coefficients  derived  from  a 
Taylor  series  expansion  of  the  forces  and  moments  acting  on  a  vessel  [5].  The  modular  modeling  approach  treats  the 
components  of  the  forces  and  moments  separately.  For  example,  the  forces  due  to  the  rudder  are  composed  of  the 
basic  forces  from  the  airfoil  in  a  flow,  the  interaction  of  the  hull  on  the  rudder  (expressed  as  flow  blockage  and 
straightening  effects),  and  the  interaction  of  the  propeller  on  the  rudder  (through  flow  acceleration  and  straightening 
due  to  the  propeller  race).  The  forces  on  the  hull  would  likewise  be  comprised  of  the  bare  hull  forces  and  the 
interactions  of  the  rudder  and  propeller  on  the  hull.  Implicit  in  the  modular  model  concept  is  the  use  of  physical 
based  models  for  each  of  the  components  [6]. 

The  modular  maneuvering  model  offers  several  advantages  over  the  coefficient  based  model.  Since  each 
component  effecting  maneuvering  performance  is  separate  in  the  modular  model,  changes  can  be  isolated.  This 
includes  changes  to  program  source  code  and  model  as  well  as  changes  to  the  ship  design.  For  example,  the  effect  of 
a  different  size  or  shape  rudder  is  isolated  to  the  rudder  model  and  the  interactions  on  other  components.  Hence,  the 
hull  force  model  or  propeller  force  model  remain  unchanged.  This  is  particularly  useful  in  preliminary  design  where 
different  concepts  are  investigated.  In  later  phases  of  design,  model  test  data  may  be  available  for  a  baseline  design. 
The  design  can  then  be  modified  with  different  propellers  or  rudders  without  necessarily  invalidating  the  previous 
model  test  work. 

By  using  the  modular  maneuvering  model  concept,  both  the  MCMSIM  and  MANSIM  program  codes  have 
been  written  in  a  highly  structured  manner  with  clear  divisions  between  force  and  moment  contributors.  This 
facilitates  better  internal  documentation  of  the  source  code  and  allows  easier  maintenance  and  future  expansion  of 
capabilities. 

This  report  documents  the  combined  maneuvering,  stationkeeping,  and  seakeeping  simulation  mathematical 
models  used  in  the  computer  program.  A  separate  programmers  manual  documents  the  computer  program  source 
code  as  implemented  for  the  MCM  and  MHC  class  vessels. 
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COMBINED  MANEUVERING,  STATIONKEEPING,  AND  SEAKEEPING 

In  developing  a  time  domain  model  of  ship  maneuvering  in  waves,  there  are  several  problems  to  be 
considered.  There  are  frequency  dependent  coefficients  in  the  equations  of  motion,  unsteady  ship  speed  and  heading, 
effects  of  the  past  history  of  a  body’s  motion  (memory  effects),  and  first  and  second  order  wave  induced  oscillatory 
exciting  forces.  Previous  attempts  to  model  ship  maneuvering  in  waves  have  generally  focused  on  the  specific 
problems  of  course  keeping  ability,  or  capsizing  and  broaching.  Few  complete  models  that  can  represent  general 
maneuvering  in  waves  have  been  developed,  with  even  fewer  implementations  of  those  models  [7].  In  general,  such 
models  are  computationally  expensive,  and  not  well  suited  to  preliminary  design  applications  or  real  time,  man-in- 
the>loop  simulations. 

The  approach  taken  here  is  to  de-couple  the  relatively  slow  maneuvering  motions  from  the  relatively  faster 
seakeeping  motions.  This  is  accomplished  by  calculating  the  calm  water  maneuvering  motions  and  adding  the  linear 
seakeeping  responses.  The  advantage  of  this  method  is  that  the  seakeeping  motions  can  be  pre-computed  outside  of 
the  maneuvering  simulation.  Figure  1  shows  the  conceptual  model  that  was  developed  as  a  framework  for  the 
development  of  the  components  of  the  model. 

It  is  important  to  note  that  when  estimating  or  experimentally  determining  the  hydrodynamic  characteristics 
of  each  component  of  the  ship  (hull,  rudders,  props,  etc.),  it  is  assumed  that  the  body  is  moving  with  sufficient  speed 
and  has  sufficient  surface  roughness  to  insure  that  the  flow  is  fully  turbulent.  However,  the  simulation  model  does 
have  the  capability  to  operate  at  very  low  speeds  where  the  flow  could  be  laminar.  In  these  cases,  appropriate 
empirical  corrections  have  been  implemented  to  provide  proper  low  speed  behavior. 

COORDINATE  SYSTEMS 

Within  this  simulator  model,  a  multitude  of  coordinate  systems  and  frames  of  reference  are  employed.  The 
standard  coordinate  systems  for  ship  motions  and  maneuvering  are  different  from  each  other  by  a  180  degree  rotation 
about  the  ship’s  x  axis.  The  bottom  of  Figure  2  shows  the  orientation  of  both  coordinate  systems  P].  In  this  work,  the 
standard  maneuvering  coordinate  system  will  be  used  unless  explicitly  stated.  Figure  2  shows  the  translating  ship 
coordinate  system  (x’,y’),  the  earth  fixed  coordinate  system  (X’,Y’),  and  the  ship  reference  coordinate  system  (x,y). 
The  positive  sense  of  the  rudder  and  drift  angles  and  propeller  and  bow  thruster  thrust  are  also  shown.  The  origin  of 
the  ship  coordinate  systems  are  at  the  ship’s  center  of  gravity.  Current  and  wind  speed  and  direction  and  predominate 
wave  direction  are  specified  in  earth  fixed  coordinate  system  using  compass  heading  angles. 
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MANEUVERING  EQUATIONS  OF  MOTION 
The  general,  nonlinear  equations  of  motion  for  a  maneuvering  ship  in  the  horizontal  plane,  with  roll,  are 
given  in  [5], 

m(u- v^)  =  ^Fx 

m(v  +  u^/)  =  ^Fy 
Izz  ij/  —  Mz 

Ixx^  +  Ixx26>„^<^  +  Ixx<»„^sin0  = 


The  cross  product  terms  on  the  left  side  of  the  x  and  y  force  equations  are  the  result  of  those  equations  being  written 
in  the  ship  coordinate  system.  Note  that  the  notation  for  the  moment  equations  has  changed  from  [1].  Rewriting 
equation  (1)  as  a  set  of  first  order  differential  equations, 

u  =  -i^(vr  +  XF>‘) 
v  =  i(-ur  +  XFy) 


r  = 


Tz7 


.  .  ,  XMx 

<!>  =  -(2<y„^<?)  +  <sin<z))  +  — 

Ixx 


(2) 


where  \j/  has  been  replaced  by  r.  The  roll  equation  is  still  shown  as  a  second  order  differential  equation  but  is 
implemented  as  a  set  of  two  coupled  first  order  differential  equations.  Relating  the  ship  coordinate  system  back  to 
the  earth  fixed  coordinate  system  yields  another  set  of  first  order  differential  equations, 

^pos  =  u  +  V  cos y/ 

Ypos  =  u  cos yz-v  siny/  (3). 

Equations  (2)  and  (3),  when  taken  together,  form  a  set  of  seven,  first  order  nonlinear  state  equations*  .  The 
general  form  is  given  by 

X  =  A(X)  +  B(t/,X)  (4) 


where  X  is  the  state  variable  vector,  U  is  the  external  force  vector,  and  A  and  B  are  matrices  representing 
the  characteristics  of  the  system.  Solution  of  the  time  evolution  of  these  equations  is  performed  using  either  a  first 
order  Euler  time  step, 


*  the  second  order  equation  for  roll  in  equation  (2)  can  be  written  as  a  set  of  two  first  order  equations. 
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(5a) 


X"*'  =X"+[A(X")  +  B(f7",X")]At 

or  a  fourth  order  Runge-Kutta  time  step, 

k,  =[A(X")  +  B(f/",Z")]At 

kj  =  [A( X "  +  (k,  /  2))  +  B(f/ " ,  Z "  +  (k,  /  2))]At 

k3  =  [A(X"  +  (kj  /  2))  +  B([/ " ,  X"  +  (kj  /  2))]At 

k4=[A(X"  +  k3)  +  B(t/",X"  +  k3)]At 

X"^'  =  X"+^[k,+2k3  +  2k3  +  kJ  (5b). 

The  first  order  Euler  method  is  generally  stable  (for  dynamically  stable  ships)  and  will  not  propagate  errors  if  the 
time  step  size  is  relatively  small  compared  to  the  speed  of  ship  maneuvering  motions.  However,  due  to  the  presence 
of  the  restoring  force  term  in  the  roll  equation,  the  first  order  integration  is  numerically  unstable  if  used.  Therefore, 
the  higher  order  method  is  required  when  roll  is  included.  If  roll  is  not  included,  the  first  order  method  may  be  used 
and  the  expense  of  the  extra  derivative  evaluations  avoided. 


EXTERNAL  FORCES 

If  we  disregard  the  seakeeping  aspects  of  the  problem  for  now,  the  total  external  forces  and  moments  acting 
on  a  maneuvering  ship  are  represented  by  ^  ^  z  ’  ^  x  •  These  summations  can  be  expanded 

into  components  for  the  hull,  rudders,  propellers,  thrusters,  wind,  and  waves  [5], 


+  Fx 

+  Fx 

+  F, 

+  F, 

+  F, 

^rudder 

props 

*  thruster 

"wind 

^  waves 

+  Fv 

+  Fv 

+  F„ 

+  F 

+  F„ 

y  rudder 

y  props 

y  Ihnistcr 

ywiid 

y  waves 

— -M  "b  M  "b^z  ”b  M  +  M 

^  ^hull  ^rudder  ^props  ^thruster  ^wind 

“  ^x  "b  M  "b  M  “b  M  +M  4"  M 

““  *  ''hull  ^^rudder  *Droos  ^thruster  *wind 


(6). 


CURRENT  EFFECTS 

Hull  forces  are  due  the  rotation  and  translation  of  the  vessel  through  the  water.  In  order  to  include  the 
effects  of  a  current,  the  relative  velocity  between  the  ship  and  the  current  must  be  calculated. 

u,  =  Vcos/?  +  V„  -  y/) 

Vc  =  -  Vsin/?  -  -  y/)  (7) 

where  the  ship  velocity  and  drift  angle  are 


5 


(8) 


V  = 

Y  pos  pos 

V 

n  “^pos 

/?=V^-atan^ 

^pos 

and  the  drift  angle  is  considered  zero  if  both  the  Xpos  and  Ypos  velocities  are  zero.  The  velocity  and  drift  of  the  ship 
relative  to  the  current  is  therefore 

5  =  atari—  (9). 

u 

Using  the  relative  velocity  through  the  water  requires  that  the  surge  and  sway  equations  in  equation  (2)  be 
rewritten  as 

K  =  m(Vcr  +  XF") 

Vc  =i(-Uer  +  XFy)  (10). 

Equation  (3)  must  then  be  corrected  for  the  current  and  written  as 

Xpos  ^  sin^  + V  cos^ 

Ypos  =  u  cos^-v  siin/  (11) 

where 

U  =  Uc-V„COS(V^-^ounem) 

V  =  V,  +  v„  sin(  y/  -  )  (12). 


HULL  FORCE  MODEL 

As  stated  in  the  introduction,  the  traditional  approach  to  modeling  hull  forces  has  been  to  assume  the 
linearity  of  forces  and  expand  them  with  a  Taylor  series  expansion.  This  is  still  the  most  straight  forward  approach  to 
developing  hull  forces  as  a  function  of  surge,  sway,  and  yaw  motions.  Here  we  will  not  include  the  effects  of  rudder, 
props,  or  bow  thrusters  in  the  expansion. 

If  the  hull  forces  are  considered  to  be  functions  of  surge,  sway,  and  yaw  motion  only  ,  we  may  write 
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^  2  <3F  2  7  ‘^x 

’‘“  d^  dw^  d[  di^  d<l>^  d<l)^^  ds  a- 

^  ^  ^  dF  <^y  2  ^y  ^y 

^yhi,  +‘3rr+-^r^  +  +^v  +  -^f+.... 

dv  Bi  di  dcf)  dcf)  dw  dt 

d!s/[  ^  2  <^Z  <^7  7  <^7  <^7  7  <^7  <^7 

“-=^''"1^'' 


.,  oM  oM  2  <^x  ^x  7  ^x  ,  <^x  7  ^x  ^x 

M,  =-^v  +  -z-^y^  +^^r+^-;rr^  +-:::^v+^r^f+. 


^hull 


dw  dw^ 


di 


d(^  ^  df 


By 


B^ 


(13) 


where  any  potential  higher  order  and  cross  coupling  terms  have  been  omitted.  The  derivatives,  also  known  as 

BF 

hydrodynamic  or  maneuvering  coefficients,  are  usually  written  in  a  shorthand  notation  such  that  — =  X„ , 


Bw 


» 


dw 


dw 


Ny  ,  and  — z~  =  ,  for  example.  The  acceleration  terms  in  the  equations  basically 

dw 


represent  low  frequency  added  mass  effects  [5]. 

Hydrodynamic  coefficients  can  be  obtained  from  captive  model  tests,  estimation,  or  system  identification 
from  full  scale  trials  or  free  running  model  experiments.  Each  of  these  methods  has  their  own  merits  and  drawbacks. 
A  discussion  of  these  is  beyond  the  scope  of  this  document  and  the  reader  is  referred  to  the  references  [5,  8,  9,  10, 
11]. 


The  normal  practice  is  to  express  the  hydrodynamic  coefficients  in  non-dimensional  form,  for  example 


(14). 


Here,  the  velocity  in  the  denominator  is  the  relative  velocity  through  the  water.  When  a  non-dimensional 
hydrodynamic  coefficient  is  used  to  calculate  the  non-dimensional  force  or  moment,  it  implies  that  the  associated 
variable  (v  in  the  example  case  of  equation  (14))  to  be  multiplied  is  also  in  its  non-dimensional  form.  Hence,  the 
dimensional  force  is  obtained  from  summing  the  products  of  the  non-dimensional  coefficients  and  associated  non- 
dimensional  variables,  and  then  multiplying  by  the  factor  ^ (L^  for  moments) 

Within  the  series  expansion  for  the  x  force,  a  term  has  been  added  to  account  for  the  hull  resistance.  The 
effective  power  (EHP)  curve  is  represented  as  a  polynomial  curve  from  which  the  hull  resistance  is  obtained, 

Pe(V)  =  V'[c,V  +  C2V'+C3V'] 


Rj 


550Pe 

V 


(15). 


The  polynomial  coefficients  are  obtained  from  curve  fits  of  upright  straight  line  (zero  drift  angle)  resistance  tests 
[12].  The  ship  speed  applied  to  equation  (15)  in  the  program  implementation  is  the  x  component  of  the  relative  speed 
through  the  water  (Uc).  In  the  special  case  of  zero  ship  speed,  the  resistance  is  zero. 
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When  running  the  simulator  program  for  low  ship  speeds,  it  was  observed  that  the  behavior  of  the  model 
was  incorrect.  If  the  vessel  was  allowed  to  decelerate  from  some  nominal  forward  speed,  the  model  did  not  show  the 
vessel  coming  to  rest  within  a  reasonable  time.  To  correct  this  behavior,  an  empirical  correction  has  been  made  for 
ship  speeds  of  5  knots  and  below.  This  correction  takes  the  form  of  a  multiplying  factor  to  the  hull  resistance  Rj  in 
equation  (15).  Table  1  lists  the  correction  factors  used  in  the  model. 

RUDDER  FORCE  MODEL 

In  the  case  of  the  MHC,  the  rudders  are  replaced  by  the  omni  directional  thrust  provided  by  the  vertical  axis 
propellers.  However,  for  the  MCM,  the  rudder  is  the  primary  device  used  to  maintain  directional  control.  In  order  to 
determine  the  forces  acting  on  the  rudder,  the  actual  velocity  and  direction  of  the  flow  at  the  rudder  must  be 
determined  from  the  general  motion  of  the  ship.  Figure  3  shows  a  schematic  representation  of  a  ship  in  a  turn.  If  the 
velocity  components  of  the  ship  center  of  gravity  (eg)  and  the  yaw  rate  are  known  then  the  drift  angle,  turning  radius, 
and  inflow  angle  at  the  rudder  (or  any  arbitrary  point  along  the  vessel  centerline)  is  calculated  from 


=  sign(r)  •  +(4g  -  x^^sinA) 

y^rud  =y^-sign(r)atan 


vLcg-XrudSinA 


Vrud  =|rL™d| 

«  =  Aud+^-'^  06). 

In  the  case  of  r=0  and  V^O,  then  the  turning  radius  is  infinite  and  (3^^  =  and  Vnid=Vc.  The  rudder  angle 
is  given  by  the  commanded  rudder  angle  up  to  the  point  that  a  set  maximum  rudder  rate  or  the  maximum  allowed 
angle  is  achieved. 

When  the  point  in  question  is  not  on  the  ship  centerline,  an  additional  correction  is  required.  If  the  turning 
radius  and  drift  angle  of  a  point  on  the  centerline  is  Lcl  and  pcL>  then  the  turning  radius  and  drift  angle  of  a  point  off 
the  centerline  is 

Lp,  =  sign(LcL)  •  ^(yptSin^L^  +  (Lcl  -  YptCOS^cL)' 


A>,  =A:L  +  sign(LcL)atan 


yp,sin/?cL 


(17). 
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The  correction,  e  in  equation  (16),  is  a  flow  straightening  coefficient  which  represents  the  blockage  of  the 
hull  to  the  flow.  It  has  been  seen  that  the  flow  straightening  is  primarily  dependent  on  the  geometric  drift  angle  at  the 
rudder,  [13].  A  relation  derived  from  the  data  presented  in  [13]  is  employed, 

f  =  Co+C,y0„,+cX+C3y5^d  (18) 

where  it  is  assumed  that  e=0  if  |  >  90  degrees. 

There  is  an  additional  interaction  effect  is  on  the  effective  flow  velocity  when  the  propeller  race  is  directed 
over  the  rudder  [  14], 

=  ^  [{>  -  »t)v....cosA^  +  kV„f  +  -  Wt)"  (19) 

where 

V„  =  -(l  -  Wt  )v„,cos^„,  +  V(l- w  +|:K.rn'D' 

k  =  0.5  +  0.26527tanh(l  .2775^)  +  0.17533tanh(2.555^) 

^=2Xp.,/D  (20). 

Some  portion  of  the  propeller  race  will  be  directed  over  the  rudder  when:  Uc>0  and  p/d>0  ,  and  possibly  Uc<0  and 
p/d>0. 

The  forces  acting  on  the  rudder  are  given  by  the  lift  and  drag  coefficients.  As  the  rudder  will  operate  in  a 
diversity  of  flow  conditions,  lift  and  drag  characteristics  are  needed  for  forward,  reverse,  and  sideways  flow  [15]. 
Like  the  hull  hydrodynamic  coefficients,  rudder  lift  and  drag  coefficients  are  generally  determined  in  a  fully 
turbulent  flow  condition.  These  coefficients  are  used  at  all  times  in  the  simulator  model  for  both  port  and  starboard 
rudders,  even  for  low  speed  operation. 

Lift  and  drag  are  defined  respectively  as  parallel  and  normal  to  the  flow.  In  order  to  obtain  the  forces  acting 
on  the  hull,  a  transformation  is  performed  from  the  flow  reference  frame  to  the  rudder  reference  frame  and  finally  to 
the  moving  ship’s  coordinate  system.  This  is  done  using 

=  CLCOsar  +  CjjSino: 

=C^cosa-C^sina  (2i) 

and 
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Fx  =(CL,,cos^-CD„^sin^)(ipV^,Aj 
Py  =(-CD.,cos^-CL^^sin<5)(i/7ViA,) 


M.=F^x„,±^F, 


M 


X 


=  FyZ„, 


(22). 


For  the  yaw  moment,  the  positive  sign  is  taken  for  the  port  rudder  and  the  negative  for  the  starboard  rudder. 
Though  the  rudders  on  the  MCM  operate  together,  each  will  experience  a  different  inflow  geometry  and  potentially 
different  propeller  races.  Therefore,  the  forces  and  moments  should  be  calculated  for  each  rudder  separately  and 
summed  together.  The  machinery  model  for  the  steering  gear  imposes  limits  on  the  extremes  of  position  and  the 
maximum  rate  of  change  of  position. 


PROPELLER  FORCE  MODEL 


The  propellers  are  the  primary  means  for  the  ship  operator  to  control  the  speed  of  his  vessel  and  for  the 
MHC,  they  are  also  the  only  devices  for  providing  directional  control.  For  the  MCM,  the  propellers  can  provide 
some  directional  control  when  the  two  propellers  are  operated  individually  at  different  speeds  and/or  pitch  settings. 

Since  the  MCM  has  controllable  reversible  pitch  (CRP)  propellers,  the  full  four  quadrant  open  water 
propeller  performance  at  all  pitch  settings  must  be  modeled.  The  open  water  propeller  data  is  represented  as  curves 
of  thrust  and  torque  coefficients  (Kj  and  Kq  respectively)  over  a  range  of  negative  and  positive  advance  coefficient, 
J  [16].  These  are  defined  as 


nD’ 


J  —  _  .  Ky  —  y  ^21^4  >  ~ 


(23) 


where  Va  is  the  speed  of  advance  through  the  water,  n  is  the  propeller  speed,  D  is  the  propeller  diameter,  T  is  the 
propeller  thrust,  and  Q  is  the  propeller  torque.  For  the  purposes  of  this  model,  the  speed  of  advance  is 

V3=V(1-W,)  (24) 


where  (1-Wt)  is  the  thrust  wake  fraction  and  V  is  the  x  component  of  the  velocity  of  a  point  at  the  propeller  hub 
(refer  to  the  discussion  in  the  previous  rudder  force  section).  The  thrust  wake  fraction  represents  a  blockage  effect  of 
the  hull  on  the  flow  to  the  propeller  and  is  usually  determined  by  powering  experiments.  In  straight  line  motion  (zero 
drift  angle),  the  wake  fraction  is  the  measured  value  from  experiments  [12].  However,  as  the  ship  maneuvers  and 
operates  at  non-zero  drift  angles,  there  will  be  less  blockage  by  the  hull  on  the  flow  to  the  propellers.  For  inflow 
angles  less  than  90  degrees,  a  cosine  squared  correction  is  used, 

(l  -  Wt)  =  +  (l  -  (25) 
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and  (1“Wt)=1  for  angles  greater  then  90  degrees.  Here  is  the  drift  angle  of  a  point  at  the  propeller  hub.  Since  the 
propellers  are  to  operated  independently  from  each  other,  the  velocity  V  in  equation  (24)  and  drift  angle  pprop  in 
equation  (25)  are  calculated  for  both  port  and  starboard  propellers. 

A  polynomial  curve  fit  for  each  Kt  and  Kq  curve  is  made  for  each  pitch/diameter  value  and  both  positive 
and  negative  J.  Reference  [17]  gives  twelfth  order  polynomial  coefficients  for  the  MCM  design  propeller.  In  order  to 
obtain  the  operating  values  for  either  port  or  starboard  propeller,  the  polynomials  are  evaluated  at  the  appropriate 
value  of  advance  coefficient.  Linear  interpolation  is  used  between  curves  to  the  actual  propeller  pitch/diameter 
setting.  The  engine  model  dynamics  and  state  equations  (discussed  later)  determine  the  rpm  to  be  used  in  calculating 
J  for  the  port  and  starboard  props. 

Having  the  total  thrust  calculated  from  the  port  and  starboard  propeller  Kj  values,  the  thrust  must  be 
decreased  by  the  thrust  deduction  factor.  Thrust  deduction  is  corrected  for  non-zero  inflow  angles  in  a  similar  manor 
to  that  used  on  the  wake  fraction  in  equation  (25).  Propeller  thrust  acts  only  along  the  longitudinal  direction  of  the 
ship  hull.  No  oblique  flow  effect  are  included  except  for  the  wake  fraction  and  thrust  deduction  factor  corrections. 

As  a  result,  and  by  assumption,  no  side  forces  are  exerted.  For  the  MCM  model,  the  forces  and  moments  delivered 
by  the  propellers  to  the  ship  are  given  by, 


F,=K^(l-tXi/7n^D'') 

F,=0 


=  0  (26). 


For  the  yaw  moment,  the  positive  sign  is  taken  for  the  port  propeller  and  the  negative  for  the  starboard  propeller. 

The  MHC  vertical  axis  propellers  use  a  different  model  based  on  information  provided  by  J.M.  Voith 
GmbH*  .  Define  the  thrust  and  torque  coefficients  (Ks  and  Kp  respectively)  as. 


*  Unpublished  data. 
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Ups  = 


'U 


ps 


Ks  = 


LbiatieU| 

2Q 

KpI^orbi.LMadeU 


K.= 


(27). 


The  value  of  Ups  is  the  orbital  velocity  of  a  propeller  blade  (tangential  velocity)  and  X  is  the  ratio  of  the  inflow 
velocity  at  the  propeller  to  the  blade  orbital  velocity.  Note  that  for  the  data  obtained,  X  is  limited  in  value  to  the 
range  of  0  to  1.  Open  water  characteristics  of  the  vertical  axis  propellers  are  given  as  curves  of  Ks  and  Kd  versus  X 
for  various  pitch  settings.  Polynomial  curve  fits  are  used  to  determine  the  values  of  Ks  and  Kd  for  given  values  of  X 
and  pitch  setting.  In  a  manner  similar  to  that  used  in  the  rudder  model,  the  force  produced  by  the  propeller  and 
imparted  to  the  hull  must  be  resolved  into  longitudinal  and  lateral  forces.  For  the  MHC,  the  forces  and  moments 
delivered  by  the  propellers  to  the  ship  are  given  by, 

«  =  y9prop+^ 

Iprop-x  ~  orbit  Lyade  Ups  )( 1  “  t)  COS  Of 

Fprop-y  =  KsCK^^orbi-UbladeU^sKl-OsinOf 
Fx  =  Fprop-x  cos^op  -  Fprop.y  siny^pr^p 
Fy  =  Fprop-y  cos^p^pp  +  Fp^pp.,  siii^p^pp 


M  =Fx  +-^^p 

IVlz  i^yApj.Qp  —  2 

M,  =  FyZp^pp  (28). 

For  the  yaw  moment,  the  positive  sign  is  taken  for  the  port  propeller  and  the  negative  for  the  starboard  propeller. 
Here  the  forces  are  first  resolved  into  xy  components  relative  to  the  inflow  velocity  and  then  resolved  into  ship  xy 
components.  Note  that  a  flow  straightening  correction  is  not  included  in  the  definition  of  a. 

BOW  THRUSTER  FORCE  MODEL 

The  bow  thruster  is  an  auxiliary  device  for  providing  directional  control.  The  MCM  class  of  vessels  use  a 
bow  thruster  which  has  a  sea  chest  in  the  keel  of  the  vessel  for  a  water  inlet  and  a  split  so  that  the  flow  travels  to 
either  a  port  or  starboard  side  outlet.  For  the  purposes  of  this  simulator,  a  relatively  simple  model  will  be  used.  Note 
that  the  MHC  does  not  have  a  thruster  of  any  kind. 
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The  ship  operator  will  issue  a  thruster  command  of  between  +100%  and  -100%  thrust.  A  second  order 
polynomial  is  used  to  transform  the  command  setting  to  the  full  scale  thrust.  The  data  needed  to  perform  the  curve  fit 
is  usually  provided  by  the  thruster  manufacture  and  should  assume  no  forward  speed.  Forward  speed  degradation 
effects  have  been  studied  by  Chislett  and  Bjorheden  [18]  and  McCreight  has  performed  curve  fits  of  the  published 
data  [14].  The  thruster  forces  can  be  determined  from 

F  =0 

X 

Fy  =  T[exp(- 1 3.3m' )  +  0.627396m  -  0.385772m'  +  0.1 24873m'  ] 


^  thruster 


M  =  F  z 

J.VX  ^  1  y  Zi 


where 


m  = 


(29) 


(30) 


Ml 


and  T  is  the  thruster  thrust  with  no  ship  velocity.  It  is  assumed  that  m=0  if  T=0.  For  a  ship  moving  aft,  it  is  assumed 
that  there  is  negligible  degradation.  This  is  represented  here  by  m=0  for  Pc=180  degrees.  The  velocity  and  drift 
angle  at  the  thruster  should  also  be  calculated  for  the  ship  in  a  turn  using  the  method  outlined  in  equation  (16). 

WIND  FORCE  MODEL 

In  a  manner  similar  to  the  current  and  relative  speed  through  the  water,  a  transform  between  the  wind  speed 
and  direction  in  the  earth  fixed  reference  system  and  the  ship  reference  system  is  given  by 

u„  =  Vcos^  +  cos(97„i„,  -  i/r) 

=  -Vsin/?-  -  y/) 


V  =  Ju^  + 

w  V  w  w 


=  atan 


-  v„ 


u. 


(31). 


The  sway  force  and  moment  coefficients  are  defined  as  the  side  force  and  yaw  moment  (non- 
dimensionalized)  acting  on  the  ship  when  the  wind  is  blowing  across  the  port  beam.  The  surge  force  and  moment 
coefficients  are  defined  when  the  wind  is  blowing  across  the  bow.  The  total  force  and  moment  contribution  from  the 
wind  is  a  blended  value  of  the  two  sets  of  wind  force  and  moment  coefficients  [14], 
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(32). 


Fx  =  iP^yl  cosy?,  ) 

Fy  =i/^airVw(AreayCy^^^sin;5,) 

M,  =  \p^yl  [(Area.C,,^^^^  cos2y?, )+ (Area^Cy,^^^^  sin2y?, )] 

M  =  F  z 

X  y  windage 


WAVE  DRIFT  FORCE  MODEL 

The  wave  drift  forces  being  addressed  here  are  the  mean,  second  order  wave  drift  forces.  This  is  an 
acceptable  model  for  use  in  real  time  man-in-the-loop  simulation  where  calculation  of  the  full  time  varying  second 
order  forces  would  be  burdensome.  The  mean  wave  drift  forces  are  represented  by 

Fx  “  PS®  Cl/3^D,  ("^0  ’  ^^wave  ) 

Fy  =  PgB  ^/3C  (Tq  ,  ) 

M  z  ~  PS®  D,  (Tq  5  ^^wave  ) 

M,  =FyZ,,,,  (33) 


where 


(34) 


is  the  non-dimensional  wave  modal  period  [14]. 

McCreight  and  Jones  have  investigated  the  wave  drift  forces  on  an  MCM  model  in  uni-directional  regular 
seas  [19].  Since  then,  McCreight  has  developed  an  unpublished  interpolation  method  for  the  original  model  test  data. 
The  non-dimensional  wave  drift  force  coefficients  used  the  simulation  computer  program  are  determined  by 
interpolation  of  the  test  data. 


ENGINE  MACHINERY  MODEL 

The  MCM-1  is  powered  by  four  Waukesha  LN  1616  DSIN  diesel  engines  geared  to  two  Transamerica- 
Delaval  reduction  gears.  There  are  also  twin  light  load  electric  motors  that  can  be  used.  For  this  model,  it  is 
reasonable  to  assume  that  two  diesel  engines  driving  one  shaft  are  represented  by  one  single  large  engine  driving  its 
own  shaft  with  no  cross  coupling  between  shafts.  No  separate  model  for  the  light  load  motors  is  included.  Even 
though  the  MHC  power  plants  are  different  than  those  on  the  MCM,  the  basic  power  plant  dynamics  represented  by 
this  model  are  sufficient  for  performing  simulations  for  both  vessels. 

The  original  engine  machinery  model  used  in  MCMSIM  Ver.  1  was  a  constant  torque  model  [20].  While 
that  model  worked,  it  did  require  solving  for  the  roots  of  a  twelfth  order  polynomial.  Experience  with  the  constant 
torque  model  implemented  in  MCMSIM  Ver.  1  has  shown  that  in  transient  conditions,  the  actual  torque  output  of  the 
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engine  (computed  from  root  solving)  exhibited  unrealistic  behavior  (signal  spikes).  This  behavior  led  to  the 
development  and  inclusion  of  a  more  realistic  and  physical  based  engine  model. 

Define  the  rate  of  change  of  the  propeller  shaft  speed  to  be, 

.  Q.-Qp^ 

- - j - 

**  shaft 


(35) 


where  Qe  is  the  torque  produced  by  the  engine  for  a  given  throttle  setting,  Qprop  is  the  torque  absorbed  by  the 
propeller,  and  Jshaft  is  the  rotational  inertia  of  the  propeller,  shaft,  ge^s,  and  engine  crank  shaft.  This  is  a  first  order 
differential  equation  and  is  included  in  the  total  system  of  equations  solved  in  equation  (4).  The  torque  absorbed  by 
the  propeller  is  computed  from, 


Qprop 


[forMCM] 

[forMHC] 


(36) 


depending  on  which  vessel  is  being  simulated. 

The  engine  throttle  setting  is  given  as  0  to  100%  where  maximum  engine  torque  is  ordered  for  a  100% 
throttle  setting.  However,  engine  torque  does  not  respond  instantaneously  to  throttle  setting.  Therefore,  a  critically 
damped  second  order  control  system  is  used  to  model  the  response  of  the  engine  to  throttle  setting.  If  the  actual 
engine  throttle  setting  is  T  and  the  commanded  throttle  setting  is  Xc,  then  the  engine  torque  is  computed  from, 

Qe  =(Zgear^shaf,)(^eMAx)  (37). 


The  second  order  control  system  is  implemented  as  a  pair  of  coupled  first  order  differential  equations  and  solved 
with  the  equations  of  motion.  Losses  in  the  system  are  accounted  for  by  the  gear  ratio  Xgear  and  the  efficiency  of  the 
system  Tjshaft  (includes  losses  due  to  friction).  For  the  throttle  control  system  to  be  critically  damped,  ^throt  =  1*  The 
system  can  then  be  defined  in  terms  of  a  “rise  time”  which  is  the  time  that  it  takes  for  the  throttle  position  to  rise 
from  10%  to  90%  setting  (starting  from  0%)  due  to  a  step  command  of  0%  to  100%  [21], 

s3.5/t,,,  (38). 


SEAKEEPING  EFFECTS 

The  six  degree  of  freedom  responses  are  governed  by  frequency  domain  transfer  functions,  calculated  using 
a  linear  ship  motions  computer  program  [22].  Time  histories  for  the  six  degrees  of  freedom  are  generated  from  the 
amplitude  and  phase  information  contained  within  the  transfer  functions, 

(t)  =  S  yos(o},^  t  +  7,  +  )]  (39). 

k 
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This  gives  the  motion  at  time  t  where  the  summation  is  over  discrete  frequency  values,  Ra  and  ep  is  the  transfer 
function  amplitude  and  phase  at  the  wave  frequency,  cOg  is  the  frequency  of  encounter  at  the  wave  frequency, 
and  is  the  wave  height  at  the  k*^  wave  frequency  defined  by 


(40). 


A  uniformly  distributed  random  phase  angle,  Yk,  is  included.  Equation  (39)  applies  specifically  to  long  crested  seas. 
Short  crested  sea  responses  require  an  additional  summation  over  wave  direction  [23]. 

The  time  history  data  files  used  are  composed  of  576  individual  six  degree  of  freedom  time  histories.  They 
are  for  six  ship  speeds  (0  to  10  knots  in  2  knot  increments),  twenty  four  headings  (0  to  345  degrees  in  15  degree 
increments),  and  four  significant  wave  height/modal  period  combinations;  1  foot  and  7  seconds  (Sea  State  2),  3  feet 
and  7  seconds  (Sea  State  3),  6.2  feet  and  9  seconds  (Sea  State  4),  and  10.7  feet  and  9  seconds  (Sea  State  5).  Each  of 
these  individual  time  histories  is  10  minutes  in  length  at  4  samples  per  seconds  with  short  crested  seas  assumed.  All 
of  the  individual  time  histories  are  concatenated  with  the  others  in  a  particular  order  to  create  a  single  composite  time 
history  data  file.  Each  composite  time  history  data  file  therefor  represents  a  10  minute  sample  of  ship  motions  for  a 
particular  sea  state. 

Access  to  the  time  history  data  by  the  maneuvering  simulator  is  performed  using  a  data  table  look-up  with 
linear  interpolation  between  heading  and  speed  as  needed.  However,  it  should  be  understood  that  linear  ship  motion 
theory  assumes  constant  mean  heading  and  constant  mean  ship  speed.  Both  the  transfer  functions  and  generated  time 
histories  are  sensitive  to  speed  and  heading  variations.  Hence,  an  interpolated  time  history  for  a  given  mode  of 
motion  may  not  necessarily  be  the  same  as  a  time  history  derived  from  a  transfer  function  at  that  identical  speed  and 
heading. 

Interpolation  between  speed  and  heading  is  performed  using  linear  blending  functions.  These  functions  are 
dependent  only  on  speed  or  heading  and  can  therefore  be  used  for  all  modes  of  motion  without  recomputing.  For 
example,  the  ship  speed  and  relative  heading  to  the  waves  are  given  at  discrete  values  Vi,...,Vn  and  ^],...,A.n 
respectively.  The  blending  functions  for  ship  speed  (<t>)  and  heading  ('F)  when  interpolating  to  speed  V  and  heading 
X  between  the  i^  and  the  i+1^  speed  and  j*  and  j+1^  heading  are  then, 

V- V 

^i+1 

a)(i)=l-<D(i  +  l) 

+  (41). 


Interpolating  at  time  t  takes  place  between  time  history  values  at  speeds  i  and  i+1  and  heading]  and  j+1, 
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(42). 


i7(t,V,.l)=[?;(e,V,,„/l,.,)4'(j  +  l)+)7(t,V„4,)>F(j)](I>(i  +  l)+ 

Kt.  V„  4,  + 1)+  r^t,  V, ,  I,  )'P(j)]4>(i) 

SUMMARY 

A  description  of  the  mathematical  model  used  in  developing  a  simulator  for  ship  maneuvering, 
stationkeeping,  and  seakeeping  has  been  presented.  The  model  is  based  on  the  concept  of  a  modular  ship 
maneuvering  model.  This  treats  the  ship  hull,  rudders,  propellers,  and  bow  thruster  as  individual  components, 
modeled  separately  using  methods  appropriate  for  each  component  with  interaction  effects  between  components 
accounted  for. 

Two  of  the  originally  planned  modifications  to  the  model  have  been  incorporated;  the  inclusion  of  ship  roll 
to  the  maneuvering  equation  of  motion  and  the  ability  to  model  vertical  axes  propellers.  Within  the  computer 
program  implementation  of  the  model,  these  two  modifications  led  to  the  complete  reformulation  of  the  structure  of 
the  generation  of  the  state  equations,  the  addition  of  a  higher  order  integration  method,  and  a  rewrite  of  the  engine 
dynamics  model.  The  implementation  of  MANSIM  Version  2  has  also  included  a  graphical  presentation  of  the 
simulation  as  it  advances  in  time. 

More  changes  to  the  computer  model  and  computer  program  code  are  planned.  These  changes  could 
include;  towing  forces,  shallow  water  effects,  ship-ship  interaction  effects,  ship-shore  interaction  effects,  time 
varying  second  order  wave  drift  forces,  man  in  the  loop  control  input,  and  automatic  control  systems  for  dynamic 
positioning.  An  important  change  for  the  computer  program  will  be  the  replacement  of  the  “hard  wired”  input  data, 
which  represent  the  vessel  model,  with  a  file  based  data  representation.  This  will  make  the  computer  program 
capable  of  more  general  vessel  simulations.  Another  important  modification  is  to  improve  upon  the  method  used  for 
the  superposition  of  seakeeping  motions  on  the  calm  water  maneuvering  motions,  possible  adding  the  capability  to 
solve  the  full  maneuvering  in  waves  problem. 
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Table  1.  Low  Speed  Hull  Resistance  Correction  Factors 


Ship  Speed  (knots) 

Hull  Resistance  Goto 

Vk  >  5.0 

1.0 

5.0  >  Vk  >  4.0 

1.2 

4.0  >Vk  >  3.0 

2.5 

3.0  >  Vk  >  2.0 

5.0 

2.0  >Vk  >  1.0 

10.0 

1.0  >Vk  >  0.5 

50.0 

Vk  <  0.5 

100.0 
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Coordinate  System  and 
Sign  Conventions 


(Earth  fixed  coordinate  system) 


Figure  2.  Coordinate  Systems  and  Sign  Conventions 
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